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To  study  the  transient  response  of  a  diffusion  flame  to  an  unsteady  flowfield,  quantitative  measurements 
of  velocity,  using  particle-imaging  veloeimetry,  and  OH  measurements,  using  planar  laser-induced  fluores¬ 
cence,  were  made  simultaneously  in  an  oscillating  counterflow  diffusion  flame.  These  non-intrusive  mea¬ 
surements  were  performed  to  spatially  and  temporally  resolve  flowfield  and  flame  characteristics  as  a 
function  of  initial  steady  strain  rate  and  forcing  frequency.  For  the  forcing  frequencies  considered  in  this 
study,  the  strain  rate  fluctuations  were  found  to  lag  the  velocity  fluctuations,  but  the  phase  difference 
decreased  with  increasing  forcing  frequency.  At  lower  forcing  frequencies,  the  width  of  the  OH  field 
responded  quasi-steadily,  but  as  the  forcing  frequency  increased,  the  OH  field  showed  transient  effects. 
The  dilatation  velocity,  defined  as  the  difference  between  the  minimum  velocity  in  the  preheat  zone  and 
the  maximum  velocity  in  the  reaction  zone,  was  used  as  a  flame  temperature  indicator.  The  dilatation 
velocity  revealed  that  the  phase  difference  between  the  velocity  and  the  temperature  increased  with  in¬ 
creasing  forcing  frequency,  confirming  the  existence  of  a  diffusion  limited  response.  The  results  presented 
here  help  to  illuminate  the  interconnecting  relationships  between  the  chemistry,  fluid  dynamics,  and  re¬ 
actant  transport  times. 


Introduction 

Turbulent  diffusion  flames  are  of  particular  inter¬ 
est  because  of  their  presence  in  most  practical  com¬ 
bustion  devices.  Flamelet  theory  is  a  method  that 
characterizes  turbulent  diffusion  flames  as  a  collec¬ 
tion  of  strained,  laminar,  one-dimensional  flamelets 
[1,2]  which  are  solely  dependent  upon  the  mixture 
fraction  and  instantaneous  scalar  dissipation  rate. 
The  flamelets  are  assumed  to  respond  quasi-steadily 
to  the  unsteady  strain  rates  of  the  turbulent  flow- 
field.  Recent  computational  and  experimental  stud¬ 
ies,  however,  have  shown  that  conditions  exist  where 
this  assumption  is  invalid  [3-6],  If  the  turbulent  Rey¬ 
nolds  number  is  sufficiently  large,  there  exists  a 
range  of  eddy  sizes  where  the  characteristic  turnover 
time  of  the  smallest  eddies  is  comparable  to  the  dif¬ 
fusion  time  of  the  laminar  flamelet  [3],  The  large- 
scale  eddies  establish  the  magnitude  of  the  mean 


strain  rate,  while  the  fluctuations  around  this  mean 
value  are  caused  by  the  smaller  eddies.  These  con¬ 
ditions  can  lead  to  a  wide  range  of  characteristic  fre¬ 
quencies  [4],  Thus,  it  is  necessary  to  investigate  the 
frequency  response  of  flamelets  in  an  attempt  to  ex¬ 
tend  the  applicability  of  the  flamelet  theory. 

The  purpose  of  this  investigation  was  to  quantify 
the  strain  rate  and  relative  OH  concentration  fluc¬ 
tuations  in  a  propane/air  flame  subjected  to  velocity 
fluctuations.  A  counterflow  diffusion  burner  was 
used  because  it  has  the  same  scalar  structure  as  a 
flamelet  in  the  mixing  zone  of  a  turbulent  reacting 
flowfield  [1].  Particle-imaging  veloeimetry  (PIV) 
measurements  were  used  to  quantify  the  velocity, 
strain  rate,  and  phase-angle  relationships  between 
measured  parameters  as  a  function  of  velocity  oscil¬ 
lation  frequency.  Simultaneously,  the  OH  field  was 
measured  using  planar  laser-induced  fluorescence 
(PLIF).  Reaction  zone  thickness  as  measured  by  the 
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LAMINAR  DIFFUSION  FLAMES— Counterflow  Flames 


TABLE  1 

Forcing  Amplitude  Shown  as  a  Fraction  of  Flow 
Reversal  and  Extinction  Voltages  from  Ref.  [8] 
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width  of  the  OH  field,  rather  than  OII  concentra¬ 
tion,  is  presented  in  this  paper. 


Experimental  Apparatus  and  Procedure 

A  schematic  of  the  counterflow  diffusion  flame 
burner  used  in  this  study  is  shown  in  Fig.  1.  The 
burner  is  a  modification  of  the  design  described  by 
Puri  and  Seshadri  [7]  and  described  in  detail  in  Ref. 
[8],  Briefly,  the  oxidizer  and  fuel  tubes  are  25.4  mm 
in  diameter  and  have  a  12.7  mm  separation  distance. 
Plenums  located  on  the  fuel  and  oxidizer  sides  are 
capped  with  20  cm  loudspeakers,  which  impose  the 
velocity  fluctuations  and  are  driven  by  an  amplified 
signal  generator. 


An  unsteady  flowfield  was  imposed  on  the  coun¬ 
terflow  diffusion  flame  by  inputting  a  sinusoidal  volt¬ 
age  signal  to  the  speakers.  The  velocity  and  relative 
[OH]  measurements  were  made  as  a  function  of  ini¬ 
tial  steady  strain  rate  (SSR)  and  forcing  frequency. 
These  measurements  were  made  at  four  temporal 
locations  within  the  sinusoidal  voltage  oscillation  ap¬ 
plied  to  the  speakers:  the  zero  amplitude  with  posi¬ 
tive  slope  (0  +  ),  maximum  amplitude  (Max),  zero 
amplitude  with  negative  slope  (0  — ),  and  minimum 
amplitude  (Min).  The  forcing  frequencies  consid¬ 
ered  in  this  study  were  30,  50,  100,  and  200  Hz.  A 
weak  forcing  amplitude  was  applied  to  the  flow;  Ta¬ 
ble  1  defines  this  amplitude  relative  to  flow  reversed 
and  extinction  as  previously  measured  [8], 

Figure  2  is  a  schematic  of  the  PIV-OH  measure¬ 
ment  system.  The  PIV  measurements  were  made 
using  two  frequency-doubled  Nd:YAG  lasers  (d  = 
532  nm);  in  the  probe  region,  the  thickness  of  the 
two  sheets  were  set  to  500  /an  with  a  height  of  20 
mm.  Hollow  zeospheres,  having  a  mean  diameter  of 
2.2  /an,  were  used  to  seed  the  airflow.  The  PIV  im¬ 
ages  were  recorded  using  a  Kodak  interline-transfer 
digital  camera  (1008  X  1012  pixels)  fitted  with  a  105 
mm  lens,  using  an/  stop  of  8.  An  interference  filter 
centered  on  532  nm  with  a  full  width  at  half  maxi¬ 
mum  (FWHM)  of  3  nm  was  used  to  block  back¬ 
ground  flame  emission.  Interrogation  regions  com¬ 
posed  of  64  pixel  squares  (0.9  X  0.9  mm)  with  75% 
overlapping  were  employed  for  the  velocity  mea¬ 
surements.  Custom  software  including  a  cross-cor¬ 
relation  algorithm  was  used  to  derive  the  velocities. 

The  OH  measurements  were  made  using  a 
Nd:YAG  pumped  dye  laser,  which  was  running  Rho- 
damine  590  dye,  and  the  output  of  which  was  fre¬ 
quency  doubled  down  to  ~282  nm,  yielding  a  10  ns 
pulse  of  approximately  10  mj  in  a  sheet  of  approxi¬ 
mately  300  /mi  in  width  and  20  mm  in  height.  The 
R/8)  transition  of  the  +  <-  X277(l,0)  band  was 
excited,  partially  saturating  the  transition,  and  de¬ 
tected  using  a  Princeton  Instruments  Intensified 
Charge-Coupled  Device  (ICCD)  camera.  The  Ri(8) 
transition  was  chosen  because  of  its  reduced  sensi¬ 
tivity  to  temperature  fluctuations;  indeed,  the  max¬ 
imum  OH  concentration  varied  by  less  than  5%  from 
the  average  cycle  value  for  any  of  the  experimental 
conditions.  A  105  mm//4.5  lens  fitted  with  UG-11 
and  WG-295  filters  was  used  to  collect  the  OII  fluo¬ 
rescence. 

The  response  of  the  reaction  zone  to  velocity  os¬ 
cillations  was  characterized  by  measuring  the 
FWHM  of  the  OH  field.  The  strain  rate  was  deter¬ 
mined  by  evaluating  the  gradient  of  the  axial  velocity 
profile  at  the  centerline  of  the  burner  on  the  oxidizer 
side  of  the  stagnation  plane.  The  centerline  velocity 
was  determined  by  averaging  the  three  closest  ve¬ 
locity  vectors  on  both  sides  of  the  centerline  that 
were  at  the  same  axial  distance  from  the  air  tube  exit. 


SIMULTANEOUS  MEASUREMENTS  IN  UNSTEADY  CFDF 


2023 


Cylindrical  Lens  &  ND  Filter 


Counterflow  Burner 

Fig.  2.  Optical  layout  for  PIV  and  OH-PLIF  measurements. 
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Fig.  3.  Velocity  field  data  showing  the  linear  fitted  line 
for  strain  rate  calculations.  Also  shown  are  the  velocities 
used  to  determine  the  dilatation  velocity. 
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This  averaged  the  velocity  over  1.16  mm  in  the  trans¬ 
verse  coordinate  at  each  axial  position,  which  was 
less  than  5%  of  the  air  tube  exit  diameter.  A  line  was 
fit  to  the  velocity  data  prior  to  the  preheat  zone,  and 
the  strain  rate  was  then  calculated  from  the  fitted 
line.  The  location  of  the  fitted  line  relative  to  the 
preheat  zone  is  illustrated  in  Fig.  3.  Twenty  OH  and 
PIV  images  were  taken  and  averaged  at  each  of  the 
four  temporal  locations  within  the  speaker  oscilla¬ 
tion.  No  corrections  for  thermophoretic  effects  have 
been  performed  on  the  velocity  measurements,  as 
previous  studies  have  shown  that  accurate  velocity 


measurements  can  be  performed  prior  to  the  pre¬ 
heat  zone  without  accounting  for  thermophoresis 


Results  and  Discussion 

Due  to  the  stoichiometry  in  a  counterflowing  pro¬ 
pane/air  diffusion  flame,  the  flame  is  located  on  the 
air  side  of  the  stagnation  plane.  The  transport  of  fuel 
is  a  result  of  diffusion,  while  the  transport  of  oxidizer 
to  the  flame  is  a  result  of  convection  and  diffusion. 
This  convective-diffusive  coupling  on  the  air  side  ef¬ 
fectively  links  the  chemistry  of  the  flame  with  the 
flowfield  dynamics.  A  flame  that  is  subject  to  a  cy¬ 
clical  velocity  field,  with  a  characteristic  cycle  time 
much  greater  than  that  of  any  of  the  relevant  trans¬ 
port  times  of  the  flame,  should  respond  in  a  quasi¬ 
steady  manner.  As  the  characteristic  cycle  time  of 
the  velocity  fluctuations  decreases  and  approaches 
the  relevant  transport  times  of  the  flame,  unsteady 
or  transient  characteristics  in  the  flame  behavior  be¬ 
gin  to  appear. 

The  measured  strain  rates  and  FWHM  of  the  OH 
field  for  an  initial  SSR  of  74  s-1  are  plotted  in  Fig. 
4  at  the  four  temporal  locations  in  the  speaker  os¬ 
cillation  for  all  four  forcing  frequencies.  At  a  forcing 
frequency  of  30  Hz,  the  instantaneous  strain  rate 
continuously  increases  from  the  temporal  locations 
of  0  +  to  0  — ,  as  seen  in  Fig.  4a.  During  this  period, 
the  OH  field  narrows  continuously.  Between  the  0  — 
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Fig.  4.  Strain  rate  and  FWHM  of  the  OH  field  results 
for  SSR  74  •  s  1  and  all  forcing  frequencies.  Error  bars  of 
one  standard  deviation  have  been  included  for  the  OH 
fields  FWHM  on  the  30  Hz  and  200  Hz  forcing  frequencies 
and  for  the  strain  rate  at  30  Hz.  Error  bars  for  the  mea¬ 
sured  strain  rate  are  excluded  on  4d  as  well  as  other  plots 
for  clarity,  but  the  magnitudes  are  very  similar  to  that 
shown  in  4a.  A  sine  wave  curve  is  shown  in  4c  fitted  to  the 
strain  rate. 


and  Min  locations,  the  instantaneous  strain  rate  de¬ 
creases  and  the  OH  layer  thickens.  It  was  also  found 
that  when  the  instantaneous  strain  rate  for  the 
forced  flames  was  above  the  steady  strain  rate,  the 
OH  field  was  thinner  than  the  steady  case.  The  con¬ 
verse  of  this  was  also  found  to  be  true.  The  strain 
rate  and  OH  field  behavior  at  a  forcing  frequency  of 
50  Hz  mimicked  the  behavior  found  at  30  Hz,  as 


seen  in  Fig.  4b.  If  the  flame  were  responding  stead¬ 
ily,  increasing  the  strain  rate  would  be  accompanied 
by  an  increase  in  the  convective  velocities  at  the 
edge  of  the  convective-diffusive  zone  of  the  flame. 
With  this  elevated  supply  of  cool  reactants,  the  flame 
would  move  closer  to  the  stagnation  plane,  due  to 
stoichiometry,  and  the  flame  would  also  cool  and  be¬ 
come  thinner.  Thinning  of  the  flame  with  increasing 
strain  rate  has  been  found  previously  in  experimen¬ 
tal  and  numerical  studies  [10,11],  Thus,  from  Fig. 
4a  and  b,  the  OII  field,  and  hence  the  flame,  is  re¬ 
sponding  to  the  instantaneous  strain  rate  in  a  quasi¬ 
steady  manner,  at  forcing  frequencies  of  30  and  50 
Hz.  Note  that  cyclic  behavior  of  the  OH  concentra¬ 
tion  has  been  found  previously  in  a  N2  diluted  meth¬ 
ane/air  flame  [6], 

As  the  forcing  frequency  was  increased  to  100  Hz, 
the  strain  rate  was  found  to  increase  continuously 
between  the  temporal  locations  of  0  +  and  0  —  and 
then  decreased  between  the  locations  of  0—  and 
Min,  as  seen  in  Fig.  4c.  This  is  similar  to  the  behavior 
found  for  lower-frequency  oscillations.  The  OH  field 
width  was  found  to  increase  between  the  temporal 
location  of  0  +  and  Max  and  then  decreased  contin¬ 
uously  between  the  Max  and  Min  locations.  Thus, 
for  the  100  Hz  case,  between  the  0  +  and  Max  lo¬ 
cations,  there  is  an  increase  in  the  instantaneous 
strain  rate  that  is  accompanied  by  an  increase  in  the 
width  of  the  OH  layer.  Also,  between  the  0  —  and 
Min  locations,  a  decrease  in  the  instantaneous  strain 
rate  is  accompanied  by  a  decreasing  width  of  the  OH 
layer.  Similar  behavior  was  also  noted  when  the  forc¬ 
ing  frequency  was  increased  to  200  Hz,  which  is 
shown  in  Fig.  4d.  At  200  Hz,  the  strain  rate  de¬ 
creased  between  0  +  and  Max  and  then  continuously 
increased  between  the  Max  and  Min  temporal  lo¬ 
cations.  The  OH  field  width  increased  between  the 
locations  of  0+  and  0—  and  then  decreased  be¬ 
tween  0  —  and  Min.  Between  the  Max  and  0  —  lo¬ 
cations,  there  is  an  increasing  strain  rate  and  an  in¬ 
creasing  OH  field  width.  Thus,  the  OH  field  does 
not  respond  quasi-steadily  to  the  instantaneous 
strain  rate  at  forcing  frequencies  above  100  Hz. 

As  illustrated  in  Fig.  4,  the  maximum  strain  rate 
occurs  at  the  0—  temporal  location  for  the  three 
lower  frequencies,  but  at  the  Min  location  for  the 
200  Hz  case.  Although  a  frequency-dependent 
phase-angle  relationship  between  the  strain  rate  and 
width  of  the  OH  field  seems  apparent,  it  cannot  be 
directly  deduced  from  Fig.  4.  This  is  due  to  a  fre¬ 
quency-dependent  phase  angle  that  exists  between 
the  speaker  motion  and  the  applied  voltage  [8].  Al¬ 
though  only  four  temporal  locations  were  measured 
for  each  flame  condition,  phase-angle  relationships 
were  still  resolvable  with  this  limited  temporal  sam¬ 
pling.  Because  the  speakers  are  oscillating  in  a  si¬ 
nusoidal  manner,  previous  studies  have  shown  that 
parameters  dependent  upon  their  motion,  such  as 
the  air  exit  velocity,  respond  similarly  [6,8,12],  The 
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Fig.  5.  Normalized  velocity  and 
strain  rate  data  with  their  fitted 
curves  for  SSR  74  s_1  and  50  Hz. 


Velocity-Strain  Rate  Phase  Relationship 


Fig.  6.  Phase  difference  relation¬ 
ships  between  the  air  velocity  and 
the  strain  rate  as  a  function  of  fre¬ 
quency  for  initial  steady  strain  rates 
of  23  s~\  44  s-1,  and  74  s~x. 


measured  results  of  the  velocity,  strain  rate,  and  the 
OH  field  were  fitted  to  an  offset  sine  function,  as 
shown  by  equation  1,  to  determine  phase-angle  re¬ 
lationships. 

V  =  V0  +  Vj  ■  sin(co  •  t  +  V2)  (1) 

Here,  V0  is  an  offset  value,  Vi  is  half  the  peak-to- 
peak  difference,  V2  is  the  phase  difference,  f  is  the 
temporal  location  in  the  oscillation,  and  w  is  equal 
to  2nf,  where/ is  the  forcing  frequency.  The  param¬ 
eters  Vo,  Vj  ,  and  V2  were  determined  from  the  mea¬ 
sured  data. 

Figure  5  shows  a  typical  plot  of  the  normalized 
measured  velocity  and  strain  rate  with  the  corre¬ 
sponding  fitted  curves  plotted  versus  time  for  a  forc¬ 
ing  frequency  of  50  Hz.  Subtracting  V0  and  then 
dividing  by  Vi  normalizes  the  measured  data  and 
fitted  curves.  Good  agreement  was  found  between 
the  measured  data  and  the  fitted  sine  functions.  The 
velocity  measurements  reported  in  Fig.  5  corre¬ 
spond  to  the  closest  point  to  the  flame  that  was  not 
affected  by  heat  release.  This  location  was  used  be¬ 
cause  the  relevant  marker  of  the  flowfield  exists  in 
the  region  where  strain  rate  is  measured  instead  of 


the  nozzle  exit  [12],  This  is  a  result  of  a  phase  dif¬ 
ference  between  the  air  exit  velocity  and  the  veloc¬ 
ities  just  prior  to  the  preheat  zone.  As  can  be  seen 
in  Fig.  5,  the  strain  rate  is  lagging  the  velocity.  Under 
isothermal  conditions  or  cyclic  oscillations  which 
have  cycle  times  much  longer  than  relevant  transport 
times  of  the  flame,  the  velocity  and  the  strain  rate 
will  be  in  phase.  Fig.  6  shows  the  phase-angle  rela¬ 
tionship  between  the  velocity  and  strain  rate  as  a 
function  of  frequency  for  SSRs  of  23  s_1,  44  s  , 
and  74  s_1.  For  all  strain  rates,  the  phase  difference 
between  the  strain  rate  and  the  velocity  is  decreasing 
as  the  forcing  frequency  is  increased  from  30  to  200 
Hz.  Egolfopoulos  and  Campbell  [4]  performed  a  de¬ 
tailed  numerical  study  of  the  frequency  response  of 
counterflowing  strained  diffusion  methane/oxygen/ 
nitrogen  flames.  For  a  flame  with  a  mean  strain  rate 
of  400  s  ,  they  found  the  strain  rate  was  in  phase 
with  the  velocity  for  forcing  frequencies  below  10 
Hz.  Above  forcing  frequencies  of  10  Hz,  a  phase 
difference  was  present  between  the  strain  rate  and 
the  air  exit  velocity,  which  reached  a  maximum  and 
then  began  to  decay  back  to  zero  as  the  forcing  fre¬ 
quency  continued  to  increased.  The  phase  differ¬ 
ence  between  the  velocity  and  the  strain  rate  was 
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Fig.  7.  Phase  difference  relation¬ 
ships  between  the  dilatation  velocity 
and  the  air  velocity  as  a  function  of 
frequency  for  initial  steady  strain 
rates  of  23  s~\  44  s~\  and  74  s_1. 


found  to  be  induced  by  a  systematic  increase  in  a  lag 
between  the  velocity  and  thermal  expansion,  which 
follows  the  phase-shifted  maximum  cyclic  tempera¬ 
ture.  The  magnitudes  of  the  oscillations  in  maximum 
temperature  were  also  found  to  decrease  as  the  forc¬ 
ing  frequency  increased.  Sensitivity  of  strain  rate  to 
thermal  expansion  has  been  observed  and  investi¬ 
gated  in  previous  studies  [10,13,14],  Experimental 
results  of  the  phase  difference  between  velocity  and 
strain  rate  presented  in  this  study  support  the  nu¬ 
merical  results  of  Egolfopoulos  and  Campbell  [4], 

Figure  6  also  illustrates  that,  generally,  the  phase 
difference  between  the  velocity  and  strain  rate  is  de¬ 
creasing  with  decreasing  initial  steady  strain  rate.  In¬ 
creasing  mean  strain  rates  result  in  increased  mean 
velocities  at  the  edge  of  the  convective-diffusive 
zone  on  the  air  side  and  thereby  reduce  the  diffusive 
time  necessaiy  for  the  oxidizer  to  be  transported  to 
the  reaction  zone.  Consequently,  higher  forcing 
frequencies  will  be  necessaiy  to  induce  a  phase 
difference.  An  analogy  has  been  presented  be¬ 
tween  the  attenuation  of  spatial  gradients  in  vis¬ 
cosity  for  Stokes  second  problem  and  the  attenu¬ 
ation  of  spatial  gradients  of  reactants  with 
increasing  forcing  frequency  [4],  With  increasing 
strain  rate,  the  frequency  at  which  phase-shifted 
temperature  (thermal  expansion)  effects  decrease 
their  influence  on  the  strain  rate  will  increase.  Ul¬ 
timately,  this  means  an  increasing  strain  rate 
should  shift  the  location  where  the  phase  differ¬ 
ence  of  the  strain  rate  becomes  non-zero,  reaches 
a  maximum,  and  decreases  back  to  zero  to  higher 
forcing  frequencies.  The  general  trends  shown  in 
Fig.  6  support  this  description. 

A  flame  temperature  marker  was  derived  to  help 
confirm  the  phase-shifting  behavior  of  the  maximum 
temperature.  Although  no  corrections  for  thermo- 
phoretic  forces  have  been  made  to  the  velocities  in 
the  heat  release  zone,  the  difference  between  the 
maximum  velocity  in  the  reaction  zone  and  the  min¬ 
imum  velocity  in  the  preheat  zone,  as  shown  in  Fig. 
3,  was  used  as  an  indicator  of  the  flame  temperature. 
The  increasing  velocities  in  the  reaction  zone  are  a 


result  of  density  changes  resulting  from  heat  release, 
so  the  general  behavior  of  the  changes  in  velocity 
should  follow  the  temperature  trend.  This  velocity 
difference  will  be  referred  to  as  the  dilatation  veloc¬ 
ity.  Fig.  7  shows  the  phase  difference  between  the 
air  velocity  and  the  dilatation  velocity  as  a  function 
of  frequency  and  initial  steady  strain  rate.  The  mag¬ 
nitude  of  the  phase  difference  was  found  to  increase 
with  increasing  forcing  frequency,  which  suggests 
there  is  a  systematic  increase  in  the  phase  difference 
between  the  velocity  and  the  temperature.  The  time 
associated  with  the  increasing  phase  difference  is 
likely  related  to  a  diffusion  time.  As  the  forcing  fre¬ 
quency  is  increased,  the  time  rate  of  change  of  re¬ 
actants  delivered  to  the  edge  of  the  convective-dif¬ 
fusive  zone  increases:  however,  a  finite  amount  of 
time  is  still  necessary  for  the  reactants  to  diffuse  to 
the  flame  front.  As  the  forcing  frequency  increases, 
this  diffusion  time  becomes  larger  relative  to  the  cy¬ 
cle  time  of  the  oscillation,  which  in  turn  shows  up  as 
an  increasing  phase  difference.  This  diffusion-lim¬ 
ited  response  has  been  observed  in  previous  numer¬ 
ical  studies  [3,4], 

Conclusions 

The  response  of  a  propane/air  counterflow  diffu¬ 
sion  flame  subjected  to  an  unsteady  flowfield  was 
investigated  to  increase  fundamental  understanding 
of  turbulent  flame  behavior  and  to  help  qualify  the 
quasi-steady  assumption  used  in  flamelet  theory.  The 
conclusions  drawn  from  this  investigation  are  as  fol¬ 
lows: 

1.  The  width  of  the  OH  field  appears  to  be  re¬ 
sponding  quasi-steadily  to  the  instantaneous 
strain  rate  for  forcing  frequencies  of  30  Hz  and 
50  Hz;  the  calculated  phase  difference  between 
the  strain  rate  and  width  of  the  OH  field  are 
183°  and  185°,  respectively.  However,  at  forcing 
frequencies  of  100  Hz  and  200  Hz,  the  phase  lag 
of  the  OH  field  thickness  grew  to  210°  and  241°, 
respectively. 
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2.  For  the  weakly  forced  flames  presented  in  this 
report,  the  sinusoidal  response  of  the  velocity, 
strain  rate,  and  width  of  the  OH  field  correlated 
well  with  the  speaker  oscillations  at  the  given 
frequency. 

3.  A  phase  difference  was  present  between  the  ve¬ 
locity  and  the  strain  rate,  which  was  found  to 
decrease  with  increasing  forcing  frequency.  This 
apparently  results  from  the  distortions  imposed 
on  the  velocity  field  due  to  a  phase-shifted  ther¬ 
mal  expansion. 

4.  The  dilatation  velocity  was  used  as  an  indicator 
of  flame  temperature.  The  phase  difference  be¬ 
tween  the  velocity  and  the  dilatation  velocity  was 
found  to  increase  with  increasing  forcing  fre¬ 
quency.  This  increasing  phase  difference  sug¬ 
gests  the  transport  of  reactants  through  the  re- 
active-diffusive  zone  of  the  flame  is  the  limiting 
step  in  the  response  of  the  flame  sheet. 
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COMMENTS 


Sebastien  Candel,  Ecole  Centrale  Paris,  France.  One 
characteristic  time  of  a  strained  diffusion  flame  is  just  the 
inverse  of  the  strain  rate.  This  could  be  used  to  define  a 
dimensionless  frequency  and  this  might  allow  collapsing 
the  data  gathered  in  these  experiments. 

Authors  Reply.  Following  the  work  of  Egolfopoulos  and 
Campbell  [4],  a  dimensionless  Stokes  parameter  may  be 
defined  as  (pf/k)  1/2,  where/ is  the  forcing  frequency  and 
k  is  the  steady  strain  rate.  If  the  change  in  thickness  of  the 
OH  zone  through  the  oscillation,  normalized  by  the  steady 
thickness,  is  plotted  against  the  log  of  this  Stokes’  param¬ 
eter,  all  the  data  from  the  four  forcing  frequencies  (30,  50, 
100,  and  200  Hz)  and  three  different  initial  strain  rates  (23, 
44,  and  74  s_1)  collapse  very  well  onto  a  single  line  with  a 
fairly  steep  negative  slope.  The  Stokes’  parameter  ranges 
from  1  to  6  in  these  experiments;  Ref.  4  shows  a  negative 
slope  when  the  normalized  maximum  temperature  is  plot¬ 
ted  over  this  range  as  well. 

• 

Jay  Jeffries,  Stanford  University,  USA.  Could  you  quan¬ 
tify  the  uncertainty  on  temperature  and  spatial  resolution? 


It  appears  you  have  drawn  conclusions  from  temperature 
variation  of  15  K  and  spatial  (OH  FWHM)  variations  of 
tens  of  micrometers. 

Authors  Reply.  Figure  4  has  been  updated  with  the  cor¬ 
rect  y- axis  scaling;  the  figure  shown  during  the  talk  and  to 
which  you  are  referring  was  off  by  a  factor  of  3.5.  Thus, 
the  actual  maximum  delta  in  OH  field  width  throughout 
the  oscillations  range  from  280  mm  at  30  Hz  to  31  mm  at 
200  Hz.  While  these  differences  are  small,  from  3%  to  30% 
of  the  steady  width  (880  mm),  they  are  real  and  very  re¬ 
peatable.  The  axial  OH  profile  was  averaged  about  the  cen¬ 
terline  in  the  radial  direction  over  3  mm  (spatial  resolution 
of  50  mm/pixel),  and  20  OH  PLIF  images  at  each  location 
within  the  oscillation  were  averaged  to  get  the  final  OH 
field  width.  The  standard  deviation  of  the  mean  among 
these  20  images  was  typically  less  than  0.25%  and  always 
less  then  1%.  Thus,  while  our  smallest  delta  in  field  width 
is  0.6  pixels  wide,  due  to  the  radial  averaging,  we  believe 
this  is  a  real  result.  The  standard  deviation  of  the  mean  in 
temperature  is  very  similar.  So,  while  the  absolute  tem¬ 
peratures  are  no  better  than  to  within  say  100  K,  the  rela¬ 
tive  temperatures  throughout  the  oscillation  are  much  bet¬ 
ter.  Further  details  regarding  the  thermometry  data  is 
forthcoming. 


